function [ejew,spectrum]=espectro(funcion,fs)
N=length(funcion);
df=fs/N;
ejew=(-fs/2:df:(N-1)*df/2);
spectrum=fft(funcion);
spectrum=abs(spectrum)/N;
spectrum=fftshift(spectrum);
